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We revisit the problem of parameter estimation of gravitational-wave chirp signals from inspi- 
ralling non-spinning compact binaries in the light of the recent extension of the post-Newtonian 
(PN) phasing formula to order (v/c) 7 beyond the leading Newtonian order. We study in detail the 
implications of higher post-Newtonian orders from 1PN up to 3.5PN in steps of 0.5PN (~ v/c), and 
examine their convergence. In both initial and advanced detectors the estimation of the chirp mass 
(M) and symmetric mass ratio (rf) improve at higher PN orders but oscillate with every half-a-PN 
order. In initial LIGO, for a IOM0-IOM0 binary at a signal-to-noise ratio (SNR) of 10, the improve- 
■ ment in the estimation of M (r\) at 3.5PN relative to 2PN is ~ 19% (52%). We compare parameter 

estimation in different detectors and assess their relative performance in two different ways: at a 
I fixed SNR, with the aim of understanding how the bandwidth improves parameter estimation, and 

for a fixed source, to gauge the importance of sensitivity. Errors in parameter estimation at a fixed 
SNR are smaller for VIRGO than for both initial and advanced LIGO. This is because of the larger 
h-j ' bandwidth over which it observes the signals. However, for sources at a fixed distance it is advanced 

LIGO that achieves the lowest errors owing to its greater sensitivity. Finally, we compute the ampli- 
• tude corrections due to the 'frequency-sweep' in the Fourier domain representation of the waveform 

within the stationary phase approximation and discuss its implication on parameter estimation. We 
, find that the amplitude corrections change the errors in M and r\ by less than 10% for initial LIGO 

. at a signal-to-noise ratio of 10. Our analysis makes explicit the significance of higher PN order 

VjQ ' modelling of the inspiralling compact binary on parameter estimation. 
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With the advent of a new generation of gravitational wave (GW) detectors such as LIGO, VIRGO, GEO and TAMA 
(3JT), 0, we are on the eve of a new era in astronomy: Gravitational Wave Astronomy (see Ref. 0,0 for recent reviews). 
The paucity of GW sources within a detectable distance, as well as the weakness of the gravitational wave signals, 
make imperative the necessity for developing optimal data analysis techniques, both for their detection and for the 
extraction of maximum information from these signals. It is for this reason that inspiralling compact binaries, which 
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I. INTRODUCTION 



can be well modelled within the general relativistic framework, have become one of the most promising candidate 
- . 1 sources for the large and medium scale gravitational wave detectors. 

An efficient data analysis scheme involves two independent aspects: first, the theoretical computation of very high 
accuracy templates and second, the design of a detection strategy adapted to the particular signal one is looking for. 
These strategies vary according to the type of signal. Gravitational waves from inspiralling binaries are transients 
lasting for a short duration in the sensitivity bandwidth of a ground-based detector. As the binary evolves the 
waveform sweeps up in frequency and amplitude, leading to a characteristic chirp signal. As the phasing of the 
waves is known very accurately, it is possible to enhance their dctcctability by using matched filtering. Bursts of 
unknown shape, as for example from a supernova, will be probed by monitoring the power excesses in the Fourier 
or time-frequency domain, but the enhancement in the visibility of the signal is not as good as when the phasing of 
the signal is known and matched filtering can be applied. In both cases, coincident observations with a network of 
detectors would assist the detection significantly, by increasing the confidence level of detection and mitigating non- 
stationarity. Continuous sinusoidal signals, as for example from a spinning neutron star, are also detected by matched 
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filtering and the signal visibility increases as the square-root of the period for which the signal is observed. Stochastic 
signals require cross-correlation of data from two or more collocated, or geographically close by, detectors. Here, the 
stochastic signal buried in one of the instruments acts as a matched filter to dig out exactly (or nearly exactly) the 
same signal in another. However, since the filter is noisy the efficiency is greatly degraded and the visibility improves 
only as the fourth-root of the duration of observation. 

As a binary inspirals adiabatically, i.e. when the inspiral time-scale is much larger than the orbital time-scale, 
it is possible to treat the problem perturbatively and expand the general relativistic equations of motion and wave 
generation as a power series in v/c, where v is the characteristic orbital velocity of the system. This post-Newtonian 
(PN) treatment has been successful in modelling the dynamics of a binary even at the late stages of inspiral and 
used in the computation of waveforms necessary for data analysis (see for a recent review) 1 . Since radiation back 
reaction causes the orbital eccentricity e to fall-off, for small e, as e oc P 19 / 18 and the orbital radius to decay much 
more slowly r oc P 2 / 3 @, the binary orbit will essentially be circular by the time the system reaches the late stages of 
the inspiral phase. Thus, in our analysis we shall restrict our attention to the case of compact binaries in quasi-circular 
orbit, i.e. circular but for the adiabatic decay of the orbit under gravitational radiation reaction. 

A. Data analysis of the chirp signal: Matched filtering 

Among the different methods suggested for the detection of chirps from inspiralling and merging binaries, matched 
filtering (also known as Weiner filtering) is the most effective technique jfl 13 • Matched filtering consists of passing 
the detector data through a linear filter, or a template, constructed from the expected signal h(t;0). Here 8 is a 
'vector' whose components are the parameters of the template. The templates h(t; 6) generally use the restricted 
waveform where for binaries in quasi-circular orbits the phase is computed at the highest PN order available, but 
the amplitude is taken to be Newtonian, involving only the dominant signal harmonic at twice the orbital frequency. 
This is different from the complete waveform, which incorporates the PN corrections to the amplitude, arising from 
the 'plus' and 'cross' GW polarizations, and hence includes the contribution from other harmonics (both higher and 
lower) besides the dominant one. Till date, for non-spinning binaries, the restricted waveform is computed to 3.5PN 
accuracy 0,0] and the complete waveform up to 2.5PN order 0,0]. The best template is probably the one which 
consists of the phasing at 3.5PN and the amplitude at 2.5PN. Presently, both the detection and parameter estimation 
problems mainly employ the restricted PN waveform although there have been some investigations on the ensuing 
improvement achieved when corrections arising from the other harmonics are incorporated by using the complete 
waveform |l3l HEj . In this paper, we confine ourselves mostly to the restricted waveform; specific amplitude 
corrections arising from the 'frequency-sweep' are considered, however, in Sec. IIVI 

In matched filtering, the unknown set of parameters characterizing the signal are measured by maximising the 
correlation of the data with a whole family of templates which correspond to different values of the parameters. The 
parameters of the template which maximises the output of a matched filter is an estimate of the true parameters. The 
parameters of a signal measured in a single experiment will be different from the actual values due to the presence 
of noise. Parameter estimation basically aims at computing the probability distribution for the measured values of 
a signal. Given a measured value from a single experiment one then uses the probability distribution function to 
compute the interval in which the true parameters of the signal lie at a specified confidence level (see Sec. ^for a 
summary of the theory of parameter estimation). In the next Section, we discuss the types of error bounds proposed 
in the literature in the context of GW data analysis. 

B. Parameter estimation of chirp signal: Different kinds of error bounds 

In parameter estimation it is of interest to obtain the distribution of the measured values and error bounds on the 
measured values of the parameters. To this end, the starting point would be to construct the Fisher information 
matrix, the inverse of which, the covariance matrix, provides an estimate of the possible errors in the measurement 
of the parameters 0. Error bounds obtained using the covariance matrix are called the Cramer- Rao bounds |l6j . 
However, for low values of the signal-to-noise ratio (SNR) the actual errors involved may be much larger than the 
errors estimated by this method. Cramer-Rao bounds fall off as the inverse of SNR, whereas the actual errors need 
not follow this behaviour. One usefulness of the Cramer-Rao bound is that, they are asymptotically valid in the limit 
of high SNR and hence provides a basis to test all other estimates. 



1 In our nomenclature, (v/c) n corresponds to ^ post-Newtonian (PN) order. Henceforth, we shall use units in which c = G = 1. 
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An alternate, and more general, way is to estimate the errors by Monte Carlo methods [171 Il8l [l9j . In this method, 
one mimics the detection problem on a computer by performing a large number of simulations corresponding to 
different realizations of the noise in each one of them. The advantage here is that, one no longer assumes a high SNR, 
which is a crucial assumption in computing the covariance matrix. In Ref. exhaustive Monte Carlo simulations 
were carried out to compute the errors in the estimation of the parameters and the covariances among them. It used 
the initial LIGO configuration and took into account only the 1PN corrections assuming, as usual, the orbit to be 
quasi-circular. It was shown that the covariance matrix grossly underestimates the errors in the estimation of the 
parameters by over a factor of two at a SNR of 10. This discrepancy disappears when the SNR is approximately 
15 for a Newtonian filter and 25 for the 1PN case. Further, the reason for the discrepancy was explained in detail 
in Ref. [2(J. Extending the Monte Carlo simulations of Ref. 18] by the inclusion of higher order terms would be 
computationally quite expensive pof. 

More rigorous bounds (Weiss- Weinstein bound and Ziv-Zakai bound) on the parameter estimation of inspiralling 
binaries are discussed in Ref. [2l| . They compare, at the Newtonian order, the results obtained by these bounds with 
the Cramer-Rao bounds and the numerical Monte Carlo results. At large SNR, they find all theoretical bounds to be 
identical and attained by Monte Carlo methods. At SNRs below 10, the Weiss- Weinstein bound and the Ziv-Zakai 
bound provide increasingly tighter lower bounds than the Cramer-Rao bound. 

C. Parameter estimation and the phasing formula: An update 

Intrinsic parameters, like masses and spins, characterising the signal can be estimated from the data collected by 
a single detector. On the other hand, the distance to the source and its position in the sky require at least three 
geographically separated detectors forming a detector network [22, H3, 13 • Cutler and Flanagan ^ have shown 
that, to a good approximation, it is sufficient to use Newtonian waveforms in these analyses. We will not, however, 
concern ourselves with the estimation of distance in the present work. 

Cutler and Flanagan [22j initiated the study of the implications of higher order phasing formula as applied to the 
parameter estimation of inspiralling binaries. They used the 1.5PN phasing formula to investigate the problem of 
parameter estimation, both for spinning and non-spinning binaries, and examined the effect of the spin-orbit parameter 
(3 (assumed constant) on the estimation of parameters. They find that parameter estimation worsens by a factor of 
about ten because of the inclusion of (3. The effect of the 2PN phasing formula was analysed independently by Poisson 
and Will |25( and Krolak, Kokkotas and Schafer 26]. In both of these works the focus was to understand the new 
spin-spin coupling term a appearing at the second PN order when the spins were aligned perpendicular to the orbital 
plane (constant (3 and a) . Compared to Ref. [26J , Ref. |25| also included the a priori information about the magnitude 
of the spin parameters, which then leads to a reduction in the rms errors in the estimation of mass parameters. It 
was shown that the effect of the inclusion of a is less drastic than (3 and that it worsens parameter estimation only by 
a factor of order unity. In a more recent work |27[ , the implications of including the spin couplings on the parameter 
estimation and the tests of alternative theories of gravity were studied using the LISA noise curve. 

D. Summary of the current work 

Starting with a brief summary of parameter estimation in Scc.[n] we discuss in Sec. IIII Al the nature of the 'chirp' 
signals from non- spinning binaries using the 3.5PN phasing formula llpfl which is now compl etely determined following 
the recent computation of the hitherto unknown parameters at 3PN |2g, |2j|, l3Ct l-'Sll |32j, l33l 1.3 11 ] . 

We study parameter estimation using three different noise curves: advanced LIGO, initial LIGO and VIRGO. Our 
choice is motivated by the fact that initial LIGO and VIRGO are the more sensitive instruments among the first 
generation of interferometric detectors with a somewhat different combination of bandwidth and sensitivity while 
advanced LIGO is prototypical of second generation instruments currently being planned. We will use the planned 
design sensitivity curves of initial LIGO and VIRGO as in Ref. |3f| and advanced LIGO 2 as in Ref. and discuss in 
Sec. IIII El the sensitivity and span of these instruments for binary coalescences. 

As mentioned earlier, Poisson and Will J25| analysed the implications of the 2PN phasing formula on parameter 
estimation of spinning binaries |36(. However, extending this to higher orders is not possible at present since spin 
effects beyond 2PN have not yet been computed. Therefore, in this work we will follow the procedure adopted in 



2 For the sake of comparison with previous work we have also carried out our study with the advanced LIGO noise curve as in Refs. I22tl25l . 
However, most of the work reported in this study uses the advanced LIGO noise curve quoted in Ref. Q. 



4 



|25| . but consider only the non-spinning case. We study in Sec. IIII HI the effect of higher order phasing terms by 
incorporating them in steps of half-a-PN order from 1PN up to 3.5PN and examine the convergence of parameter 
estimation with PN orders. We compare the errors for the different noise curves and assess their relative performance 
in two different ways: at a fixed signal-to-noise ratio (Sec. IIIi"H)l . with the aim of understanding how the sensitivity 
bandwidth improves parameter estimation, and for a fixed source (Sec. IIII D|) . to gauge the relative importance of 
sensitivity and bandwidth. We have examined the correlation of parameter estimation results to the number of useful 
cycles |23 and the sensitivity bandwidth fSec lIIIEll . which together can explain the performance of different detectors 
with regard to parameter estimation. 

In Sec. lIVI we study the effect of the amplitude terms arising from the 'frequency-sweep' dF/dt within the stationary 
phase approximation [23 . These corrections cause the SNR (which is related to the total energy emitted by the system) 
of a given binary to vary as we go from lower to higher PN orders. The results are compared against the standard 
restricted waveform approach and should be viewed as a prelude, albeit inconsistent, to parameter estimation using 
the complete waveform. We conclude in Sec. IVl with a summary of our results, their regime of validity, limitations 
and future directions. 

Our main conclusion is that the 3.5PN phasing formula leads to an improved estimate of the binary parameters. 
For instance, in the case of black hole binaries, at a SNR of 10, the estimate of chirp mass (symmetric mass ratio), 
more specifically ln.M (hi 77), improves while using the 3.5PN phasing formula as compared to the 2PN by about 
19% (52%). Improvements are seen in all cases but are relatively smaller for lighter binaries. At a fixed SNR, 
VIRGO provides a better estimate of the parameters compared to both initial and advanced LIGO configurations 
owing to its better sensitivity bandwidth. This is true over the entire mass range and even for lower mass binaries 
for which VIRGO accumulates fewer number of useful cycles. For a fixed source, however, advanced LIGO measures 
the parameters most accurately, as expected, with VIRGO doing better than initial LIGO. Our investigation of the 
amplitude corrections from 'frequency-sweep' within the stationary phase approximation finds that the percentage 
change induced by this effect in parameter estimation is less than 10% for initial LIGO at a SNR of 10. 



II. A BRIEF SUMMARY OF PARAMETER ESTIMATION THEORY 



A firm statistical foundation to the theory of gravitational wave data analysis was laid down by the works of 
e.g. Finn and Chernoff [2UE3 and Cutler and Flanagan [52]. This Section briefly outlines the problem of parameter 
estimation relevant to this paper. Notation and treatment of this Section essentially follow Ref. |2^l25l l41 | (see also 
[t! l42l I43L Elf for further details). We restrict our discussion to measurements made by a single detector. 



A. Matched filtering 



The output of a gravitational wave detector contains both the signal and noise and is schematically represented as 

x(t) = h(t) + n(t) , (2.1) 

where x(t) is the signal registered and n(t) is the noise, which is assumed to be a stationary Gaussian random variable, 
with zero mean, i.e., 



n{t) = 0. 



(2.2) 



Here an overbar denotes the ensemble average (over many realisations of the noise or, equivalently, over an ensemble 
of detectors). Let q(t) define a linear filter and c(t) its correlation with the detector output x(t) 



c(t) 



dt' x(t')q(t + t') . 



Define a new quantity er[<7](£), such that c(t) is normalized by the square root of its variance, 



a[q](t) = 



cit) 



c 2 (t)-c(t) 



1/2 



2^J c °dfx(f)q*(f)e 2 ^ t 

ur d f s h (f)ww] 1/2 



(2.3) 



(2.4) 



where x(f) and <?(/) are the Fourier transforms of x(t) and q(t), respectively, Sh{f) is the real, one-sided power 
spectral density defined only for positive frequencies by 



1 



n(f)n*(f') = ^(f-f)S h (f), 



(2.5) 
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and n(f) is the Fourier transform of n{t) defined as n(f) = dtn(i)e 2m i t . The filtered SNR is defined by the 
ensemble average 

2$if™dfh(f)q*(f)e 2 ^ t 

[J~dfS h (f)\q(f)\' 



p[q](t) = am) = r ^ w : 1/2 . (2.6) 



An optimal filter is the one which maximises the SNR at a particular instant, say t = 0, and is given by the matched 
filtering theorem as 

where 7 is an arbitrary real constant. Thus, the SNR corresponding to the optimal filter is given by 



B. Parameter estimation 



Though we may have a prior knowledge of the form of the signal we will not know what its parameters are. 
Indeed, the parameters are to be measured in the process of matched filtering. This is achieved by maximising the 
correlation in Eq. i|2.4|) with a whole family of templates corresponding to different values of the signal parameters. 
The parameters of the filter which maximise the correlation are the measured values attributed by the analyst to the 
signal presumed to be buried in the data. These parameters need not agree, in general, with the actual parameters of 
the signal since the measured values depend on a particular realization of the detector noise. 

For a given incident gravitational wave, different realizations of the noise will give rise to somewhat different best-fit 
parameters. However, if the SNR is high enough, the best-fit parameters will have a Gaussian distribution centered 
around the actual values. 

Let 9 a denote the 'true values' of the parameters and let 9 a + A9 a be the best-fit parameters in the presence of some 
realization of the noise. Then for large SNR, errors in the estimation of parameters A9 a obey a Gaussian probability 
distribution of the form [39I ] 

p(Af)= p (°' e -5 r - A9lAr , (2.9) 

where is a normalization constant. In the above expression r a & = (h a \hb), where h a = dh/d9 a , is the Fisher 
information matrix evaluated at the measured value 9 of the parameters 9. Here, ( | ) denotes the noise weighted inner 
product. Given any two functions g and h their inner product is defined as: 

(9\h) = 2j o df m . (2.10) 

Using the definition of the inner product one can re-express r a ^ more explicitly as 

k(f)h b (f) + haif)K(f) 



ab 



Sh(f) 



df. (2.11) 



The variance-covariance matrix, or simply the covariance matrix, defined as the inverse of the Fisher information 
matrix, is given by 

E afc = (A9 a A9 b ) = (r- 1 )" 6 , (2.12) 

where (•) denotes an average over the probability distribution function in Eq. I|2.9|) . The root-mean-square error a a 
in the estimation of the parameters 9 a is 

aa A 2\l/2 



a a = {(A9 a ) 2 Y / =VS^, (2.13) 
and the correlation coefficient c ab between parameters 9 a and 9 b is defined as 

cab = (A9 a A6 b ) = ^ 
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(There is no summation over repeated indices in Eqs. I|2.13[l and (|2.14|) .) As a consequence of their definition the 
correlation coefficients must lie in the range [—1,1]. When the correlation coefficient between two parameters is close 
to 1 (or —1), it indicates that the two parameters are perfectly correlated (respectively, anti-correlated) (and therefore 
redundant) and a value close to indicates that the two parameters are uncorrelated; covariance close to 1 (or — 1) 
among parameters causes a large dispersion in their measurement. 

In our analysis we will apply the method outlined above to three prototypical systems normally considered in 
gravitational wave studies related to ground-based detectors. These include a binary neutron star system (NS-NS), a 
neutron star-black hole system (NS-BH) and a binary black hole system (BH-BH). Throughout our analysis we shall 
assume that the mass of a NS is 1.4M Q and that of a BH is 1OM . 



III. PARAMETER ESTIMATION USING THE 3.5PN PHASING FORMULA 



Having outlined the essential results from the theory of parameter estimation, we proceed to address the question of 
extracting the parameters from the chirp signal using the 3.5PN phasing formula. Our computation parallels the one 
by Poisson and Will [25| except that we confine our attention to the case of non-spinning binaries whereas Ref. [2{| 
dealt with spinning binaries. 



A. Fourier transform of chirp at 3.5PN order 

To compute the Fisher information matrix we would need the Fourier transform h(f) of the signal h(t). (Note that 
here and below / is the Fourier transform variable which should not be confused with F, the instantaneous frequency 
of emitted radiation.) Following earlier works, we employ the stationary phase approximation (SPA) to evaluate 
the Fourier amplitude of the waveform. Given a function B(t) = 2 A(t) cos 4>(t), where d \nA/dt <C d(f>(t)/dt and 
\d 2 (f>/dt 2 \ <C (d<fi/dt) 2 , the SPA provides the following estimate of the Fourier transform B(f): 

B (f) ~ _^L e *[*/(*/)-f], />o, (3.1a) 

where = 2nft - <j>(t) , (3.1b) 

and — = 2irF(t). (3.1c) 
dt 

In this equation tf is defined as the time at which F(tf) = f and ^/(i/) is the value of ^/(t) at t = it. Starting 
from the 3.5PN phasing formula in [IJ}, the Fourier transform has been explicitly calculated in Refs. j35a,^]|. This 
Fourier domain waveform, which forms the basis of our further calculations, is given by 

h(f) =.4.r 7/6 e W) , (3.2) 

where A oc A4 5 / 6 Q(angles)/.D, and to 3.5PN order the phase of the Fourier domain waveform is given by 

o N 

4 128r7u a 

' k=0 

where v = (ttM/) 1 / 3 , M = mi+m2 is the total mass of the binary, r\ = m\mil M 2 is the dimensionless mass ratio and 
D the distance to the binary. We shall find it useful in our study to deal with the chirp mass defined by M. = rf/ b M 
rather than the total mass M. The coefficients a^.'s, k = 0, . . . , N, (with N = 7 at 3.5PN order) in the Fourier phase 



7 



are given by 
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(3.4b) 

(3.4c) 

(3.4d) 

(3.4e) 

(3.4f) 



(3.4g) 
(3.4h) 



Among the coefficients above, 0:5 can be simplified further. This interesting possibility arises because of the can- 
cellation of v 5 of the 2.5PN term with that of the overall factor in the denominator of Eq. I|3.3|) . Consequently, all 
but the Inv terms in 05 are constants and can be absorbed in a redefinition of the phase 3 . Indeed, we find that all 
our estimations, except A(f> c , remain unchanged irrespective of whether we choose a§ as above, or a simplified one 
retaining only the In v term. 

In the 3PN phasing, until recently there were two undetermined parameters, A and 9, arising from the incompleteness 
of the Hadamard self-field regularisation at 3PN 4 . By dimensional regularisation A and 9 have been now determined 
in Ref. |H,|2^| and [^3, 001 El H3 respectively, completing the general relativistic compact inspiral phasing to 3.5PN 
order: A = — 3§fjj — —0.6451 and 9 = — q^q 1 — —1.28. A has also been determined by alternative approach |34j . 



Following earlier works, we choose the set of independent parameters 9 describing the GW signal to be 

6= (ln^,/ t c ,(/) c ,ln7W,ln?7), 



(3.5) 



where t c refers to the coalescence time, <j> c refers to the phase at coalescence instant, /o is a scaling frequency related 
to the power spectral density (PSD) of the detectors (see next Subsection). Note that A is taken to be one of the 
independent parameters. Computing the Fisher information matrix T a b, whose elements are given by (h a \hb) (where 
a and b are indices which run over the parameters), is the first step towards our goal. The upper cut-off in computing 
the integrals in Eq. (|2.8|l and H2.11JI is taken to be the GW frequency at the last stable circular orbit (LSO) given, for 
a test mass in a Schwarzschild spacetime of mass M, to be 



^upper -^ls 



(6 3 / 2 ^l/) \ 

We take the lower limit in the integrals to be the seismic cut-off frequency f s of the detector. 



(3.6) 



B. Sensitivity and span of LIGO and VIRGO 

We compute the covariance matrix for three noise curves to understand the effect of detector characteristics on 
parameter estimation. The noise curves used are advanced LIGO as in and initial LIGO and VIRGO as in [35| . 
We have fitted the following expression to the noise PSD of advanced LIGO given in 



Sh(f) 



So 

00, f < f s , 



5x 



111(1 



+ x 4 /2) 



(l + x 2 /2) 



f>fs 



(3.7a) 
(3.7b) 



3 We thank Luc Blanchet for pointing this out to us. 

4 The ambiguity parameter occurring at 3PN should not be confused with the set of parameters 9 a describing the GW. 
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FIG. 1: Amplitude spectrum (left panel) of initial LIGO, VIRGO and advanced LIGO together with the luminosity distance 
(right panel) at which RMS-oriented binaries would produce a SNR of 5. 



where x = f / fa, fo = 215 Hz (a scaling frequency chosen for convenience), f s = 20 Hz is the lower cutoff frequency 
[defined such that for NS-NS binaries the gain in SNR by reducing the lower limit of the integral in Eq. Ij2.8[l below 
f s is less than 1%], and Sq = 10 _49 Hz _1 . Note that the above PSD is significantly different from the advanced 
LIGO noise curve used in earlier studies. Indeed, authors of Ref. 0, |22, Hjj \^ use the PSD of advanced LIGO 
to be S h (f) = So [x- 4 + 2 + 2a; 2 ] , / > / s , and S h (f) = oo, / < / s , with x = f/f , fo = 70 Hz, f s = 10 Hz and 
So = 6 x 10 _49 Hz _1 , which has a significantly better low-frequency sensitivity than what is currently believed to 
be possible for the next generation of LIGO. Hence, we have chosen to work with the more recent estimate given in 
Eq. 

The initial LIGO noise curve from Ref. j3^ is given by 



S h (f) = So [(4.49a;)- 56 + 0.16x~ 4 - 52 + 0.52 + 0.32x 2 ] , f > f s 

= oo, / < / s , 



(3.8a) 
(3.8b) 



where again x = f/f , withf = 150 Hz, f s = 40 Hz and S = 9 x 10~ 46 Hz -1 . Finally, for the VIRGO detector the 
expected PSD is given by |3q : 



S h (f) = So [(6.23a;r 5 

= OO, / < f s , 



2x~ 



1 + X 2 } , / > fs 



(3.9a) 
(3.9b) 



where fo = 500 Hz, f s = 20 Hz and So = 3.24 x 10~ 46 Hz _1 . The amplitude spectra [i.e. the square-root of the power 
spectral densities given in Eqs. H3.7fl - H3.9|) ] °f t ne various detectors are plotted in the left hand panel of Fig. ^ 

The SNR achieved by these detectors for binaries of different masses not only depends on the distance D at which 
the source is located but also on the orientation of the orbital plane with respect to the line-of-sight. In order not to 
be biased one can consider binaries of root-mean-square (RMS) orientation and compute the SNR they would produce 
in a given detector. One can turn around the question and ask the distance at which sources of RMS orientation 
would produce a specified SNR. Indeed, the distance D at which a binary of RMS orientation achieves a SNR po is 
given by [37| 



D(M,tj) 



1 



2/ .3 



2 77 M 5 / 3 
15 



/iso(Af) f-7/3 



Sh(f) 



df 



1/2 



(3.10) 



As is well known the SNR depends only on the chirp mass M. = rff^M and not on the masses of the two bodies 
separately. The SNR is maximum for equal mass binaries (for which 77 = 1/4) and is smaller by a factor \/4r/ for 
systems of the same total mass but consisting of stars of unequal masses. The right-hand panel of Fig. 2] plots the 
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TABLE I: Convergence of measurement errors from 1PN to 3.5PN at a SNR of 10 for the three prototypical binary systems: 
NS-NS, NS-BH and BH-BH using the phasing formula, in steps of 0.5PN. For each of the three detector noise curves the table 
presents At c (in msec), A(j> c (in radians), AM/M and Ar\/r\. 





NS-NS 


NS-BH 


BH-BH 


PN Order 


At c A<f) c AM/M At)/t] 


At c A(f> c AM/M A-q/r) 


At c A(f> c AM/M Ari/rj 



Advanced LIGO 



1PN 

1.5PN 

2PN 

2.5PN 

3PN 

3.5PN 



0.3977 0.9256 0.0267% 4.656% 

0.4668 1.474 0.0142% 1.638% 

0.4623 1.392 0.0143% 1.764% 

0.5090 1.359 0.0134% 1.334% 

0.4938 1.331 0.0135% 1.348% 

0.5193 1.279 0.0133% 1.319% 



0.5959 1.261 0.1420% 7.059% 
0.7394 2.091 0.0763% 2.316% 
0.7208 1.848 0.0773% 2.669% 
0.9000 1.219 0.0686% 1.515% 
0.8087 1.131 0.0698% 1.571% 
0.9966 0.9268 0.0679% 1.457% 



1.162 1.974 1.041% 59.88% 
1.441 3.188 0.6115% 9.609% 
1.404 2.850 0.6240% 10.79% 
1.819 1.574 0.5300% 5.934% 
1.544 1.580 0.5466% 6.347% 
2.078 1.161 0.5241% 5.739% 



Initial LIGO 



1PN 

1.5PN 

2PN 

2.5 

3PN 

3.5PN 



0.3598 1.238 0.0771% 9.792% 

0.4154 1.942 0.0419% 2.768% 

0.4109 1.816 0.0423% 3.007% 

0.4605 1.650 0.0384% 2.129% 

0.4402 1.618 0.0389% 2.170% 

0.4754 1.517 0.0383% 2.099% 



0.9550 2.510 0.5217% 20.06% 

1.182 4.135 0.2850% 5.410% 

1.148 3.597 0.2903% 6.316% 

1.467 1.975 0.2491% 3.305% 

1.286 1.798 0.2554% 3.474% 

1.666 1.324 0.2456% 3.151% 



2.406 5.038 4.750% 216.2% 

2.986 8.143 2.781% 28.81% 

2.900 7.179 2.851% 32.82% 

3.836 3.119 2.351% 16.48% 

3.159 3.123 2.446% 17.94% 

4.512 1.912 2.314% 15.77% 



VIRGO 



1PN 

1.5PN 

2PN 

2.5PN 

3PN 

3.5PN 



0.1363 0.5134 0.0183% 3.044% 
0.1578 0.7981 0.0098% 1.004% 
0.1562 0.7515 0.0098% 1.085% 
0.1743 0.7045 0.0091% 0.7957% 
0.1671 0.6920 0.0092% 0.8083% 
0.1797 0.6562 0.0091% 0.7858% 



0.4906 1.069 0.1134% 5.782% 
0.6069 1.763 0.0603% 1.923% 
0.5918 1.561 0.0611% 2.215% 
0.7384 1.039 0.0541% 1.263% 
0.6632 0.9672 0.0551% 1.309% 
0.8183 0.7968 0.0536% 1.215% 



1.621 1.854 0.8745% 52.12% 

1.430 2.972 0.5095% 8.586% 

1.395 2.667 0.5199% 9.625% 

1.787 1.545 0.4417% 5.370% 

1.532 1.547 0.4552% 5.724% 

2.024 1.173 0.4369% 5.201% 



luminosity distance at which binaries of RMS orientation and consisting of stars of equal masses would produce a 
SNR of po = 5. After computing the covariance matrix we shall use this plot to study how parameter estimation 
varies in different interferometers for sources at a fixed distance. 



C. Parameter estimation using 3.5PN phasing — Fixed SNR 



In this Section, we examine how the addition of higher order terms in the phasing formula affects the parameter 
estimation of the binary. We start from the 1PN phasing formula and add terms in steps of half-a-PN order up to 
3.5PN, which is the most accurate expression currently available. We are interested in the case of non-spinning binaries 
(ignoring spin and orbital angular momentum) and hence estimate only five parameters (In A, fot c -, <t> c , ln.M, In rj). We 
calculate the elements of T a b by explicitly computing the derivatives of the Fourier domain waveform with respect to 
(w.r.t) different parameters and taking their noise-weighted inner products. The derivatives and the Fisher matrices 
are too lengthy to be displayed here. We note that Ti a = 6i a p 2 , which renders the Fisher information matrix in block 
diagonal form. Since In A is now entirely uncorrelated with all other parameters, we only consider the Fisher matrix 
calculated from the partial derivatives of h(f) with respect to the four parameters (fot c , 4> c , In M., In 77). T\i can be 
thought of as an independent block, and further calculations involving A become trivial. Finally, by inverting the 
Fisher information matrix one constructs the covariance matrix. 

First, we computed the covariance matrix using the advanced LIGO noise PSD as defined in Ref. j^], which 
facilitates a comparison of our results with those discussed in the literature. Indeed, at 1.5 PN order we found our 
results in perfect agreement with the numbers given in Table I of Ref. [22] and at 2PN order our calculation reproduces 
the results in Table V of Ref. j2|J. In both of these papers, ln/i, where p is the reduced mass, is chosen to be the 
independent parameter instead of In 77. However, the errors in these quantities are simply related by — = f so 
that the comparison is straightforward. In the rest of this paper we study only the most recent advanced LIGO noise 
PSD together with initial LIGO and VIRGO. 



10 



Sources at fixed SNR Sources at fixed distance 




10° 10 1 10° 10 1 

Mass of the binary (M ) Mass of the binary (M Q ) 

FIG. 2: Comparison of errors in the estimation of t c , M and rj for sources with a fixed SNR of 10 (left panels) with those for 
systems at a fixed distance of 300 Mpc (right panels). 

Next, let us consider the covariance matrix computed using the noise PSDs of advanced and initial LIGO, and 
VIRGO, as given in Eq. (|3.7[) - H3.9[) . The errors in the measurement of the various parameters are tabulated in 
Table U for all the interferometers and for three prototypical binaries (NS-NS, NS-BH and BH-BH), assuming a 
fixed SNR of 10 in each case. Although the SNR is fixed, different detectors might accumulate the SNR over different 
bandwidths, causing the errors to be greater or smaller compared to one another. In agreement with what one expects 
intuitively based on the bandwidth of the various detectors (cf. Fig. ^ left panel), we find the errors in the various 
parameters to be the smallest for VIRGO, followed by a factor of roughly 10-70% larger errors in advanced LIGO 
compared to VIRGO, and a factor of 3 larger errors in initial LIGO compared to advanced LIGO. 

In going from lower to higher post-Newtonian order, we find that there is an 'oscillation' of the errors in the chirp 
mass and reduced mass. However, the errors at 3.5PN are always smaller than at 2PN. The opposite oscillation is 
observed for the errors in t c : the error in t c at 3.5PN is always higher than at 2PN. The fact that the reduced mass 
and chirp mass show the same trend is due to the correlation coefficient cmjj (listed in Table ITTf) all being close to 1. 

The oscillation in the variances with PN order can be partially understood by an examination of the correlation 
coefficients between t c , M. and rj. In Table ILTl we have listed the correlation coefficients together with the errors in the 
estimation of parameters in the case of advanced LIGO for a NS-BH system for all PN orders starting from Newtonian 
but let us first discuss the trend at orders beyond the 1PN correction. From this Table we see that the estimation of M 
and 77 improves (degrades) depending on whether the correlation coefficients cm^ decrease (respectively, increase) with 
varying PN order. Similarly, the estimation of t c improves (degrades) depending on whether the correlation coefficients 
c t c M ( or 7 equivalently, c tcn ) decrease (respectively, increase) with PN order. We have also checked that the estimation 
of <t> c becomes better (worse) with PN order with reduction (respectively, enhancement) in the correlation coefficients 
c <p c M ( or c <t>c.v)- same trend is seen for other systems and detector configurations, though we do not list those 
numbers to avoid proliferation of details. The behaviour of the errors at 0PN and 1PN is not in agreement with 
this general trend because at 0PN we have only three parameters - t c , 4> c and A4. As we go from 1PN to 1.5PN 
the ambiguity function greatly changes its orientation because of the change in sign in the PN series at 1.5PN [cf. 
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TABLE II: PN variation in parameter estimation and the associated correlation coefficients for the NS-BH system for the 
advanced LIGO noise curve. 



PN Order 




Ct c rj 




Af c (ms) 


AM/M (%) 


An/n (%) 

II 1 \ ' u / 


OPN 


-0.6451 






0.2775 


0.0255 




1PN 


0.8166 


-0.8810 


-0.9859 


0.5959 


0.1420 


7.059 


1.5PN 


0.7983 


0.9280 


0.9444 


0.7394 


0.0763 


2.316 


2PN 


0.7947 


0.9239 


0.9460 


0.7208 


0.0773 


2.669 


2.5PN 


0.8145 


0.9519 


0.9309 


0.9000 


0.0686 


1.515 


3PN 


0.8001 


0.9405 


0.9333 


0.8087 


0.0698 


1.571 


3.5PN 


0.8274 


0.9608 


0.9294 


0.9966 


0.0679 


1.456 



Eq. lEHc|) and Eq. fijty ]. 

Though the PN variation of parameter estimation accuracy seems to be dominantly explained by the variation of 
the correlation coefficients, it should be borne in mind that the variances in a particular parameter is a combination of 
the covariances and the availability of a greater structure or variety in the waveforms not fully assessed in this paper. 
This will be the subject of a study we shall take up in the near future; it is important to understand in more detail 
why the errors in t c worsen at higher PN orders as it has implications in the determination of the direction to the 
source. Table III summarizes the results of this Section. It provides the percentage decrease in the errors due to the 
greater accuracy (3.5PN as opposed to 2PN) in the phasing of the waves: the reduction is the highest for a BH-BH 
binary for which the improvement in the estimation of r\ is 52% and that of M. is 19% at an SNR of 10 for the initial 
LIGO noise curve. 



D. Parameter estimation using 3.5PN phasing — Fixed source 

The focus of this Section is to understand the effect of detector sensitivity (as opposed to bandwidth) on parameter 
estimation. The results of the previous Section, wherein the errors are quoted at a fixed SNR, cannot be used to gauge 
the performance of different detectors: a more sensitive detector has a larger SNR for a given source and therefore a 
better estimation of parameters. Hence, we translate the results for the errors in parameter estimation for different 
detectors but normalized to a fixed distance instead of a fixed SNR. Since the errors associated with the parameter 
estimation are inversely related to SNR (<r cx 1 / p) , given the error cto corresponding to a known SNR po (results for 
po = 10 are quoted in Table QJ, one can calculate the error a at another SNR p (corresponding to a fixed distance, 
say, 300 Mpc) by a simple rescaling of the results listed earlier. Indeed, a = pq(Jq/p, which can be recast in terms of 
the distance to the source, using Eq. (|3.10|) . as 



a(D L ) = Pa a Q 7r 2 ^D L 



' df 



15 Jf, S h (f) 



(3.H) 



Fig. summarises the results shown in Tabled (3. 5PN entries) over the entire parameter space of interest for sources 
with a fixed SNR of 10 (left panels) and also the consequent results from the scaling in Eq. (|3.11l) for sources at a fixed 
distance of 300 Mpc (right panels). The advantage of having a greater bandwidth is revealed by looking at panels 
on the left which shows the errors in VIRGO to be the smallest, followed by advanced and initial LIGO instruments. 
Although the signal-to- noise ratios in the case of VIRGO are similar to those of initial LIGO (cf. Fig. 2] right panel), 
Fig. |21 reveals that VIRGO measures the parameters more accurately. Indeed, the errors in VIRGO are smaller than 
in initial LIGO by a factor of 2 to 4 and this is entirely as a result of VIRGO's larger bandwidth. Unlike the case of 
fixed SNR, detector performance is explicit in the plots for sources at a fixed distance. It is evident that the errors 
reduce by about 30-60 times in advanced LIGO as compared to initial LIGO. Advanced LIGO gains a factor of 10-15 
in SNR relative to initial LIGO and this accounts for most of the improvement in its parameter estimation. However, 
it also gains another factor of 3 to 4 because of its greater bandwidth. From the foregoing discussion we conclude that 
as far as parameter estimation is concerned VIRGO performs better than initial LIGO and that advanced LIGO can 
measure the parameters significantly better than what one might conclude based on the improvement over VIRGO 
in its visibility of the signals. 

A final comment: The plots on the right-hand panel of Fig. [2] are somewhat flattened as compared to those on 
the left-hand panel due to the fact that errors for sources at a fixed distance are (anti) correlated with the variation 
of SNR with mass. In other words, there are two competing effects on parameter estimation as the mass of the 



12 



TABLE III: Percentage change of parameter estimation accuracy at SNR p — 10 for non-spinning compact binaries due to im- 
proved phasing accuracy from 2PN to 3.5PN. Percentage change for the parameter a n is taken to be = 100 X (l — 4' 5PN /c r f N ) • 
Negative values imply worsened parameter estimation in going from 2PN to 3.5PN. 



Interferometer 


NS-NS 


NS-BH 


BH-BH 


t c 4> c \nM inn 


t c <frc InM In 77 


t c (j> c \nM In 77 


Adv. LIGO 
Ini. LIGO 
VIRGO 


-12.33 8.118 6.993 25.23 
-15.70 16.47 9.456 30.20 
-15.05 12.68 7.143 27.58 


-38.26 49.85 12.16 45.41 
-45.12 63.19 15.40 50.11 
-38.27 48.96 12.28 45.15 


-48.01 59.26 16.01 46.81 
-55.59 73.37 18.84 51.95 
-45.09 56.02 15.97 45.96 



binary is increased. On the one hand, estimation becomes worse since the signal spends smaller amount of time in 
the detector band and the number of cycles available to discriminate different signals goes down. On the other hand, 
as we increase the mass of the binary the SNR increases thereby aiding in discriminating between different systems. 
These competing trends cause the error in the estimation of the time-of-coalescence and symmetric mass ratio to show 
a minimum for a binary of total mass M ~ 10M Q . No such minimum is seen, however, in the case of the chirp mass. 
This is because the error in the chirp mass rises more steeply with mass than the SNR can cause it to dip. 



E. Parameter estimation and Number of useful cycles 



To investigate further the correlation of parameter estimation performance with detector characteristics we consider 
the total number of cycles in the detector bandwidth and more importantly the number of useful cycles for a particular 
detector for the three systems under consideration. The total number of cycles N total, is defined as 



-Wtotal 



dF 



1 d<f> 
2^ dF 



(3.12) 



where <fi is the phase of the GW, .Fbegin and fend correspond to the upper and lower cut-off frequencies for the 
astrophysical system under consideration. Since the phasing of the waves is a post-Newtonian expansion in the 
parameter v the total number of cycles depends on the post-Newtonian order. At the dominant Newtonian order, 
assuming that the lower frequency cutoff of the detector is much smaller compared to the last stable orbit frequency 
of the system, the total number of cycles for a binary of total mass M and mass ratio rj is given by 



N, 



total 



-5/3 



327T77 



(3.13) 



The total number of cycles goes inversely as the mass ratio being the smallest (for a given total mass) for equal mass 
binaries and is quite a sharp function of the total mass. It has an artificiality to it in that it depends on the chosen 
lower-frequency cutoff, increasing quite rapidly as the the cutoff is lowered. Moreover, iVtotai nas no information 
about detector characteristics. Motivated by these facts Ref. [§3 proposed that the detector performance can be 
better understood using the idea of the number of useful cycles N usc f u \ defined as 



N, 



useful 



™%Hf)N{f) 



7 M/) 



(3.14) 



where N(F) is the instantaneous number of cycles (i.e., the number of cycles spent at the instantaneous frequency F) 
and w(f) is the weighting function that depends on the effective noise of the interferometer and the amplitude of the 
source defined as 



N(F) 



dF/dt 



if) 



a 2 (f) 



(3.15) 



with a(f) being the 'bare amplitude' appearing in the Fourier domain waveform within the SPA, — a(f)/ V F 

and h 2 n = / Sh(f). Unlike the total number of cycles, the number of useful cycles contains information about both the 
detector and the source: it is weighted by the noise PSD of the instrument and amplitude of the source. Moreover, 
while the total number of cycles depends critically on the choice of the lower-cutoff, the number of useful cycles is a 
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FIG. 3: Left hand panel is the plot of the derivative dA r use f u i/d(ln /) against the frequency (in arbitrary normalization) for the 
three detectors. Similarly, right panel gives the number of useful cycles as a function of the total mass of the binary for the 
three detectors. 

robust estimator and it is pretty much independent of the cutoffs chosen as long as the frequency range covers the 
sensitivity bandwidth of the instrument. 

At Newtonian order, the instantaneous number of cycles is given by N(f) = 5(7rAf/)- 5 / 3 /(96?7), which clearly 
exhibits the well-known fact that irrespective of the mass of the system it is best to design a detector with a good 
sensitivity at as low a frequency as possible. The instantaneous number of cycles decreases rapidly with frequency, 
but most of the contribution to the integral in Eq. I|3.14[) comes from the region of the band where weighting function 
w(f) — a 2 {f)/h? s (f) — J 1 / 3 / Sh(f), has a minimum. As shown in Fig. [3] (right-hand panel) for binaries whose total 
mass is larger than 11M© the number of useful cycles is larger in VIRGO than the other two instruments, while 
just the opposite is true for systems whose mass is smaller than 11 M . The reason for this behaviour can be seen 
by inspecting the left-hand panel of Fig. |3] where we have plotted the integrand dN use f u i/dlog f of the number of 
useful cycles [cf. Eq. (|3~TIf) ]. A binary of total mass 100M© has its last stable orbit at F\ so ~ 43 (M/100M©) -1 Hz, 
and increases in inverse proportion to the mass for systems with lower masses. Since the integral in Ea. l3.14|) is 
terminated at -Fisoj from Fig. [3] we see that as the upper limit of the integral increases (equivalently, the mass of 
the binary decreases) at first the number of useful cycles for VIRGO begins to increase. This feature explains why 
VIRGO has more number of cycles than the LIGO instruments for binaries with greater masses. However, owing to 
their relatively narrower bandwidth (as compared to VIRGO) both the LIGO instruments quickly catch up and for 
F\ so ij 300Hz, (equivalently, a total mass of M ^ 14 M©), they have greater number of useful cycles than VIRGO. 
Thus, the relatively broader bandwidth of VIRGO is responsible for the smaller number of useful cycles at lower 
masses. 

In general, one can correlate the larger errors associated with the estimation of parameters of massive systems 
with the smaller number of useful cycles for these systems (see Table ITVj) . It may be recalled that Ref. [13 showed 
that the number of useful cycles is a good quantifier of detector performance with regard to detection issues such 
as cffectualness. However, the efficiency in parameter estimation is a combination of bandwidth and the number of 
useful cycles and not the latter alone. Thus, though VIRGO has a smaller number of useful cycles than the two 
LIGO detectors for the NS-NS system, its parameter estimation at a fixed SNR is far better because of its broader 
bandwidth. 

Following Ref. [2S|. where the effects induced in parameter estimation due to the inclusion of the 2PN term was 
understood in terms of the additional total number of GW cycles accumulated at that order, we also use a very 
similar idea to understand the PN variations in parameter estimation of Table [J But unlike |25| , we use the number 
of useful cycles instead of the total number of cycles. From Table wherein we have given the errors in chirp 
mass and symmetric mass ratio together with the contributions to the useful GW cycles from each PN order term 
in phasing, it is obvious that, in general, when the number of cycles increases in g oing from one order to another, 
errors decrease (and vice versa) suggesting a possible correlation. Further, following 25], we tested this argument by 
artificially flipping the sign of each PN order term in the phasing (keeping all lower order terms with the correct sign) 
and comparing the errors. If such a correlation exists, one would expect the trend to be reversed, as the additional 
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TABLE IV: Number of useful cycles (and total number of cycles in brackets) for different systems and different detectors 
computed using 3.5PN phasing. To compute the total number of cycles the lower cut-off is chosen to be the seismic cutoff 
frequency of each detector and the upper cutoff is the frequency corresponding to the LSO. 



Detector 


NS-NS 


NS-BH 


BH-BH 


Adv. LIGO 


284 (5136) 


60 (1111) 


14 (184) 


Ini. LIGO 


251 (1615) 


59 (330) 


12 (52) 


VIRGO 


140 (5136) 


64 (1111) 


18 (184) 



number of useful cycles accumulated reverses its sign. Indeed Case B of the table does show the opposite trend 
confirming this correlation. There is an important exception to this correlation while going from Newtonian to 1PN, 
where though the number of useful cycles increase, the parameter estimation worsens. A little thought reveals that 
another more dominant aspect comes into play at this order due to the inclusion of the new parameter rj which could 
increase the errors associated with the original set of parameters. This is confirmed by looking at the parameter 
estimation of the Newtonian and 1PN orders using a smaller set of four parameters i.e. {In A, t c , <j> c , In A4}, excluding 
In ry. We find that the percentage error in chirp mass decreases from 0.0126 to 0.0120 for NS-NS case and 0.4833 
to 0.4183 in the BH-BH case in step with the increase in number of useful cycles. However, the reason behind the 
anomalous behaviour in going from 1 to 1.5PN and 3 to 3.5PN - where despite the decrease in the number of useful 
cycles, the parameter estimation improves - is not clear from the present analysis. Thus the previous considerations 
are not sufficient to completely understand the variation of parameter estimation with the PN order. 

Based on the understanding obtained in the previous paragraph, we conclude the Section with the following com- 
ment: At present we do not have a detailed understanding of the reason underlying the variation in parameter 
estimation with PN orders since the inclusion of higher PN terms could lead to one or more of the following: (a) in- 
troduction of a new parameter (e.g. rj in going from 0PN to 1PN) leading to an increase in the variance of the existing 
parameters, (b) increase in the 'variety' of waveforms leading to a reduction in the variance, and (c) change in the 
covariance among the various parameters. Though by a critical examination of the results summarised in Tables ^ 
and ITT1 some of these effects can be seen in action, it is not easy to disentangle these individual effects and present a 
consistent quantitative picture. This, we leave to a future study. 



IV. BEYOND THE RESTRICTED WAVEFORM: AMPLITUDE CORRECTIONS DUE TO 
FREQUENCY-SWEEP AND ITS IMPLICATIONS 

In the foregoing Sections we worked with the restricted PN approximation. In this approximation the GW phase 
is taken to as high a PN accuracy as available while the amplitude is assumed to be Newtonian. Indeed, all har- 
monics, except the dominant one at twice the orbital frequency, are neglected. From Eq. (|3.1() . one can see that the 
Fourier-domain amplitude is determined by the product of the time-domain amplitude A(F) oc F 2 ^ 3 and the factor 
(dF/dt)^ 1 / 2 , where dF/dt is the 'frequency-sweep' or 'chirp rate' of the signal. The frequency-sweep provides a way 
of (partially) computing the dependence of the wave amplitude on different PN orders. This correction, in addition 
to being calculable, should be of some relevance when we compared in Sec. Ill Til parameter estimation accuracy at 
different PN orders where, following Ref. [25|. we assumed the SNR to be the same at all PN orders. Our assumption 
was justified since in the restricted PN approximation there is no change in the amplitude of the signal as we go 
from one PN order to the next. However, the frequency-sweep causes the Fourier amplitude to change across the 
PN orders and leads to variations in the SNR with the PN orders. Since the errors depend on the SNR, one should 
rescale the errors by the ratio of SNRs to compare fairly the PN trends in parameter estimation of the chirp signal. 
In what follows, we will set up the necessary formulas to normalize the errors to the same SNR. However, it is imme- 
diately obvious that a more consistent calculation should begin with the full amplitude corrections arising from the 
GW polarizations computed in Ref. [TTL IT^. in lieu of the restricted approximation used here, and by including the 
sub-dominant harmonics. Inclusion of these terms is beyond the scope of this paper and will be addressed elsewhere. 

To estimate the amplitude corrections due to the frequency-sweep F, we start from the Fourier domain waveform 
in the stationary phase approximation which can be written as 



/oc 
h(t) e - 27,lft dt = 
-00 



M 

2r? — Q (angles) 



F(v) 



e lW \ (4.1) 



where v — (nMf) 1 / 3 . Using the expression for F at the Newtonian order, it can easily be shown that Eq. 14. 1|) 
reduces to Eq. I|3.2|l . From Eq. 14.1fl it is clear that the PN corrections in the frequency-sweep F [see Eq. (|4.8|) below] 
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TABLE V: Correlation of parameter estimation and number of useful cycles with PN order (n) for NS-NS and BH-BH Binaries 
for initial LIGO noise curve. Case A corresponds to the standard PN coefficients in the phasing formula (e a = 1, a < n). Case 
B refers to the results corresponding to a flip in sign of the a — nPN term keeping all other lower orders with correct sign 
(e a = — 1 for a = n and 1 for a < n). Errors listed are all in percentages. The values for the Newtonian order are obtained 
using a set of four parameters, {In A, t c , 4>c, In M}, excluding lnr;. 





NS-NS 


BH-BH 


PN Order (n) 


AM/M Ar)/r) N usefu i 


AM/M Arj/ V N usefu i 


Case A 


OPN 


0.0126 




247.8 


0.4833 




14.98 


1PN 


0.0771 


9.792 


27.13 


4.750 


216.2 


7.283 


1.5PN 


0.0419 


2.768 


-22.98 


2.781 


28.81 


-9.148 


2PN 


0.0423 


3.007 


-1.197 


2.851 


32.82 


-0.496 


2.5PN 


0.0384 


2.129 


2.406 


2.351 


16.48 


1.850 


3PN 


0.0389 


2.170 


-1.735 


2.446 


17.94 


-1.971 


3.5PN 


0.0383 


2.098 


-0.151 


2.313 


15.75 


-0.236 


Case B 


1PN 


0.0771 


8.858 


-21.65 


4.750 


158.2 


-3.579 


1.5PN 


0.0547 


1.842 


80.65 


3.237 


28.04 


21.56 


2PN 


0.0415 


2.564 


-53.64 


2.727 


25.63 


-19.30 


2.5PN 


0.0515 


5.085 


-4.700 


14.96 


473.7 


-2.395 


3PN 


0.0380 


2.089 


6.563 


2.271 


15.26 


6.461 


3.5PN 


0.0395 


2.248 


-3.453 


2.625 


20.89 


-4.978 



introduces a related PN correction in the amplitude as discussed earlier in the Section. To proceed further we note 
that the formula for F can be normalized w.r.t. its Newtonian value Fj^ and written as the product of the Newtonian 
value and PN corrections Fq: 



F = F M F C . 



(4.2) 



Schematically Fq can be written as 
F c = 



1 + F C 1PJV + F}- 5PN + F? PN + F?- 5PN + F 3PN + F?- 5PN 



Using Fa/ 



!)(> 



(nAiF) 11 / 3 and v = (nMF) 1 ^ 3 and some simple algebra, one can write, 
hc{f) = B M B c e^\ B M = Ar V \ Be = -L 



(4.3) 



(4.4) 



where Byv, as m Eq. I|2.8|l . is the Newtonian functional dependence. Using Eq. i|4.4|). the expression for SNR can be 
re- written as 



p 2 = 4 / df B c 



2 

Sh(f) 



(4.5) 



From the definition of SNR, Eq. (|2.8|) . it is clear that the SNR varies with the PN order of Fq. Similarly, one can 
write down the components of the Fisher matrix T a b as 



r„f, — 2 



df 



Bm 2 

s h (J) 



dBc dBc „ 2 di> 9V 

+ &C TT^ 



d0 a d9 b 



d0 a d9 b 



(4.6) 



where 8 a and 9 b are the parameters in the GW signal. (Be is a PN series in v and its 9 dependence arises solely from 
the mass dependence of v). In Sec lIIII Be was effectively taken to be unity. Here we relax that assumption by taking 
into account the PN corrections involved. 



16 



TABLE VI: Variation in SNR with the PN order due to the amplitude corrections arising from the frequency-sweep dF/dt in 
the stationary phase approximation. We assume initial LIGO noise spectral density and place the source such that at OPN 
order we have a SNR of 10. 



PN Order 


NS-NS 


NS-BH 


BH-BH 


OPN 


10.00 


10.00 


10.00 


1PN 


10.53 


11.26 


12.19 


1.5PN 


10.08 


9.560 


9.357 


2PN 


10.06 


9.483 


9.178 


2.5PN 


10.11 


9.736 


9.861 


3PN 


10.08 


9.355 


9.157 


3.5PN 


10.07 


9.341 


9.060 



The frequency-sweep appearing in Eq. 1)4.1(1 above can be straightforwardly calculated from the expressions for the 
flux function T and the energy function E, determining the GW phasing in the adiabatic approximation. It is given 
by M 



F(v) 



3v 2 F(v) 
'irM 2 E'iv) 1 



(4.7) 



where v — {ttM f) 1 / 3, and E' = 4-^. Using the 3.5PN accurate expression for E and T available in [T(j, the expression 
for F up to 3.5PN is given by 



dF 
dt 



3.5PJV 



90 



■(nMF) 



11/3 



5irM 2 

( 34103 13001 



59 



V 18144 
'10447322203 10tt 2 



2010 ' 18 ' 1 



4159tt 189tt 



139708800 



541 



5005 



850 



890 ' 2592 ' 105 
4415 358075 



1712 

7 - 

105 ' 

log (16 x) 



072 8 
273811877 451tt 2 



Tj ) (irMF) 5/3 
„ 010 



(ttMF) 



1088040 

2 



48 



9 



A ry 



4032 0048 



91495 2 \ . 
' 7+ 1512 77 )<« MF ) ' 



(4.8) 



•) and the coefficients A = — ill? ~ —0.6451, 6 = 



In the above expression 7 is the Euler's constant (7 = 0.577 • 

11831 no 
9240 — p_f 

In Table ED we summarize how the SNR varies with the PN order for different sources assuming that the SNR 
corresponding to the Newtonian order is 10. The convergence of the SNR's with PN orders is pretty obvious, although 
it should be recalled that the complete waveform includes PN corrections from other harmonics that are comparable 
to the higher order terms in the frequency-sweep [Til Il2j | . It would be interesting to see how the results change when 
these are included. We also note that the variation of the SNR is greater for systems with larger masses. Using 
a 3.5PN frequency-sweep, instead of the Newtonian one, increases the SNR by 0.7% for a NS-NS binary, while the 
SNR decreases by 9.5% for a BH-BH binary. Though these amplitude corrections may not be important for NS-NS 
binaries, they might be relevant for the BH-BH case. 

Using the results in Table IvTI one can implement a simple procedure to obtain better error estimates. One can 
scale the results of Sec. IIIII obtained within the restricted waveform approximation, by the factor p n /po, where p n 
and po are the SNRs at nPN and OPN orders, respectively. In this simple estimate one is effectively neglecting the 
contributions to the Fisher matrix from the variation of the F terms in the amplitude w.r.t the signal parameters 
(see Eq. ((4.0(1 . We incorporate this contribution in a more general and rigorous way in what follows. 

Our more general procedure is based on Eqs. 1(4.5(1 and 1(4. 6|) which accounts for the SNR and the Fisher matrix, 
respectively, with the full F dependence in amplitude. The steps leading to the final results listed in Table IVTT1 are 
as follows: (i) compute the amplitude A such that the SNR at OPN is 10; (ii) compute the Fisher matrix taking 
into account the amplitude corrections from the frequency-sweep using Eq. 1(4.6(1 : (iii) scale the final results by p n / Po- 
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TABLE VII: Parameter estimation with amplitude corrections from the frequency-sweep incorporated for the initial LIGO noise 
curve and SNR p = 10. nPN refers to the choice of n th PN order both in the amplitude and phase of the frequency-domain 
waveform. 



PN order 



1PN 

1.5PN 

2PN 

2.5PN 

3PN 

3.5PN 



NS-NS 



At c A(f> c AM/M Arj/v 



0.3208 1.157 0.0752% 9.431% 
0.4417 2.018 0.0426% 2.851% 
0.4427 1.903 0.0432% 3.115% 
0.4693 1.672 0.0387% 2.149% 
0.4768 1.709 0.0398% 2.261% 
0.5186 1.615 0.0393% 2.199% 



NS-BH 



At c A<j) c AM/M Arj/ri 



0.9077 2.449 0.5241% 19.94% 
1.234 4.232 0.2840% 5.484% 
1.206 3.695 0.2892% 6.408% 
1.496 2.009 0.2480% 3.326% 
1.387 1.919 0.2556% 3.583% 
1.787 1.435 0.2459% 3.260% 



BH-BH 



At c A<f>c AM/M A-q/T) 



1.717 3.544 

2.929 7.883 

2.830 6.897 

3.801 3.093 

3.106 3.068 

4.474 1.939 



3.400% 152.6% 
2.647% 27.71% 
2.680% 31.24% 
2.316% 16.28% 
2.307% 17.20% 
2.203% 15.33% 



The covariance matrix obtained from such a procedure can then be compared with that obtained in Sec. IIIII The 
procedure above is obviously equivalent to choosing a 'running' amplitude A n , with po — 10. 

In Table IVlTl the variation of errors with different PN orders is shown for the initial LIGO noise curve 5 . The 
oscillation of errors with PN orders remains after the inclusion of the frequency-sweep and one infers that changes 
due to these amplitude terms are not very significant. At an SNR of 10 the difference is at most 10%. 



V. CONCLUSION 



A. Summary and discussion of results 



We have carried out a detailed study to understand the implication of 3.5PN phasing formula on parameter es- 
timation of non-spinning binaries using the covariance matrix. We also compare parameter estimation using three 
different noise curves, advanced LIGO, initial LIGO, and VIRGO. The results of our study can be summarised as 
follows: 

1. The parameter estimation of non-spinning binaries improves significantly, as expected, by employing the 3.5PN 
phasing formula instead of the 2PN one. It is no surprise that the same trend is observed for all the three 
detectors. Improvements are larger for NS-BH and BH-BH systems and least for NS-NS binary. For initial 
LIGO, at a SNR of 10, the improvement in the estimation of parameters M. and r\ for BH-BH binaries is as 
large as 19% and 52%, respectively, whereas for NS-BH binaries it is 15% and 50%. Improvements in the case 
of VIRGO are slightly less compared to LIGO (cf. Table Hill). 

2. In proceeding from 1PN to 3.5PN, one sees an oscillation of variances with each half PN order. However, 
the errors in the mass parameters at 3.5PN are always smaller than at 1PN and one can see a convergence 
within this limited sequence. The oscillation of errors is a characteristic feature of the PN approximation. In 
Ref. [4?|. a similar oscillatory behaviour is seen in the context of the detection problem. The variation in 
parameter estimation accuracies with PN orders seem to be dominantly determined by the covariances between 
the parameters t c , cf> c , M. and r\. 

3. For sources at a fixed distance the errors in the estimation of parameters are least for advanced LIGO and the 
highest for initial LIGO, the performance of VIRGO being in between. Although initial LIGO and VIRGO 
obtain similar SNRs for sources with the total mass in the range [1, 50] M Q , the errors in VIRGO are smaller 
than in initial LIGO by a factor 2-4 due entirely to its greater bandwidth of observation. 

4. The number of useful cycles is greater in VIRGO than LIGO for higher mass binaries (M ~ 10M©) but the 
opposite is true for lower mass binaries. 

5. Parameter estimation is better if the number of useful cycles is higher but the performance also depends on the 
sensitivity bandwidth of the instrument. The notion of number of useful cycles together with bandwidth can be 
used to gauge detector performance with regard to parameter estimation. 



5 The numbers listed in Tables IVII and IVHI are those obtained by numerically integrating Eqs 14.51 and 14.61 without any further re- 
expansion of Be in Eq. 14.41 . 
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6. The variation of the Fourier amplitude of the gravitational waveform across different PN orders arising from its 
dependence on the frequency-sweep dF/dt, and its implication on parameter estimation is examined. We present 
a Table showing how the SNR varies across the PN orders for the initial LIGO noise curve. This correction 
affects the errors associated with parameter estimation by less than 10% and motivates an analysis using the 
complete waveform including all other harmonic contributions to the GW amplitude from the 'plus' and 'cross' 
polarisations which are now available up to 2.5PN in the comparable mass case |T^ . 



B. Limitations, Caveats and Future directions 



We conclude by pointing out the regime of validity of our analysis of error bounds, its limitations and possible 
future directions. 

1. Our estimates are based on the Cramer- Rao bound which is valid only in the regime of high SNR. Though at a 
SNR of 10 our calculations may be reasonably secure, in general they are less rigorous and provide only an upper 
bound on the errors involved. A full-fledged Monte-Carlo simulation would provide tighter bounds, though that 
would be computationally quite expensive. 

2. In Sec II VI we addressed the effect of inclusion of amplitude corrections arising from the frequency-sweep. This 
treatment is not fully consistent as one neglects the amplitude corrections from the other harmonics of the orbital 
frequency; a future study should address this issue more consistently. 

3. Based on the recent runs of the GW detectors LIGO and VIRGO more 'realistic' noise curves are now available. 
The parameter estimation using these realistic noise curves should be eventually addressed. 

4. A similar study in the case of spinning binaries is not possible until the terms corresponding to the effect of spins 
in the phasing formula are available beyond the present 2PN accuracy. 

5. A more detailed study is needed for completely understanding of the reasons for PN variations of the errors. We 
leave this for future study. 

6. The higher order phasing terms could play a major role also in the estimation of distance of the binary for a 
network of detectors. We will address this problem in a future work. 

While finalising this paper we learnt that E. Berti and A. Buonanno have also looked at the estimation of parameters 
using the 3.5PN phasing formula [Zsj . 
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